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Abstract: We present a calculation of heavy-flavor production in hadronic collisions in asso- 
ciation with one jet matched to parton shower Monte Carlo programs at next-to-leading order in 
perturbative QCD. Top-quark decays are included and spin correlations in the decay products are 
taken into account. The calculation builds on existing results for the radiative corrections to heavy- 
quark plus one jet production and uses the POWHEG BOX for the interface to the parton shower pro- 
grams PYTHIA or HERWIG. A broad phenomenological study for the Large Hadron Collider and the 
Tevatron is presented. In particular we study — as one important sample application — the impact of 
the parton shower on the top-quark charge asymmetry. 
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1. Introduction 

The Large Hadron Collider (LHC) and the Tevatron provide a unique experimental environment for 
top-quark physics. Akeady with the currently accumulated luminosity at both colliders, but even 
more so with the anticipated high statistics data in the LHC 7 TeV run, the prospects for precision 
studies of top-quarks are excellent. Precise experimental measurements for top-quark production 
demand theoretical predictions with comparable precision. Within Quantum Chromodynamics 
(QCD) this requires the knowledge of the hard scattering process beyond the leading order (LO) 
in perturbation theory. Predictions for the hadroproduction of top-quark pairs, i.e., pp tt, are 
available at next-to-leading order (NLO) since long [1^], and can be improved further thanks 
to resummation or the use of approximate next-to-next-to-leading order (NNLO) results, see e.g., 
[5, 6]. The NLO QCD corrections for the processes pp tt+jet [7, 8] and pp tt + ljets [9, 10] 
are more recent achievements. 

For the direct comparison with experimental data, however, e.g., in order to model experi- 
mental acceptances, it is important to provide also theory predictions for the production of fully 
exclusive events. Due to the inherent limitations of the fixed-order approaches, in the past these 
were simulated by shower Monte Carlo (SMC) programs, which were accurate to LO only. Nowa- 
days, much better theoretical precision can be reached by merging NLO computations with parton 
showers. The merging method has been pioneered by MC@NLO [11]. As an attractive feature the 
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approach combines the accuracy of exact hard matrix elements for the large angle scattering in- 
cluding the radiative corrections to first order in the strong coupling constant tts with the soft and 
collinear emission described by the parton shower. The former ingredient usually displays a re- 
duced sensitivity with respect to variations of the renormalization and factorization scales /Ir and 
/ip, while the latter accounts for the correct Sudakov suppression of colUnear and soft emissions. 
The matching between the regions of hard and of soft and collinear emissions is smooth. 

The follow-up development has seen an extension of the merging procedure by means of a 
method, dubbed POWHEG [12], which is independent of the parton SMC generator used. More- 
over, the POWHEG method allows for the generation of positive weighted events only and thus for 
a very efficient event generation. The POWHEG BOX [13] now provides a general framework for 
implementing NLO calculations in SMC programs (see also [14] for progress towards further au- 
tomation in MC@NLO). Processes involving the hadronic production of top-quarks have been subject 
to these improvements from very early on. Top-quark pair-production and single-top production 
processes have been implemented in MC@NLO [15, 16] and POWHEG [17-19] since quite some time 
and are now widely used by experimental collaborations. 

In this article, we are concerned with tt pair-production in association with one hard jet, be- 
cause a large fraction of inclusive ?f-events does actually contain one or even more additional jets. 
Moreover, due to the larger phase space, the relative importance of data samples with ff+jets in- 
creases at the LHC with respect to the Tevatron. Top-quark pair-production associated with jets 
is also an important background to Higgs boson production in vector boson fusion and for many 
signals of new physics, e.g., those motivated by supersymmetry. For the process tt+']Qi we are 
merging the NLO QCD corrections [7, 8] (see [20] for the impact of top-quark decays) with the 
parton SMC HERWIG [21] and PYTHIA [22,23] using the POWHEG BOX. An independent implemen- 
tation of this process using the results of HELAC-NLO [9,24] has been reported recently in ref. [25]. 
In this work we also include spin-correlation effects in the produced events. 

The outline is a follows. In Sec. 2 we discuss the implementation details in the POWHEG BOX, 
which are specific to the process tt+iel. We briefly describe the vaUdation procedure and care- 
fully examine details of the event generation in POWHEG. In Sec. 3 we present a phenomenological 
analysis for the Tevatron with y/s = 1.96 TeV center-of-mass energy (Sec. 3.1) and the LHC for 
^/s = 7 TeV (Sec. 3.2). In case of the Tevatron we investigate the impact of the parton shower on 
the forward-backward charge asymmetry which is currently measured at the Tevatron for the inclu- 
sive sample. For the LHC the impact on the charge asymmetry visible in the rapidity distribution is 
analyzed. Furthermore we discuss the inclusion of the top-quark decay and its impact on differen- 
tial distributions of the decay products. Our conclusions and our outlook on future developments 
are given in Sec. 4. 

2. Implementation 

The POWHEG BOX [13] provides a well-defined framework for implementing general NLO calcula- 
tions in parton SMC programs. To that end, it requires as input particular information about the 
individual components of the NLO calculation under consideration, i.e., about the Bom process, its 
virtual radiative corrections and the real emission contributions. Combination of the latter two and 
cancellation of the emerging soft and collinear singularities also requires the definition of a subtrac- 
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tion scheme. Within the POWHEG BOX the so-called Frixione-Kunszt-Signer (FKS) subtraction [26] 
is the method of choice. 

For the reaction tt+]et we are building on the computation of [7, 8] for the QCD corrections to 
NLO, which provides us with all necessary details. The Ust of all flavor structures of the partonic 
Bom and the real emission processes can be obtained from generic matrix elements by considering 
all possible crossings of Ught particles into the initial state. In Bom approximation these are given 
by 

O^ttggg, O^ttqqg, (2.1) 

and for the real corrections by 

O^ttgggg, O^ttqqgg, O^ttqqq'q', O^ttqqqq, (2.2) 

with q,q' = {u,d,s,c,b} and q / q'. 

We use MadGraph [27] for the computation of all Born squared amplitudes from eq. (2.1) and 
real emission matrix elements squared from eq. (2.2). Likewise, the color correlated Bom squared 
amplitudes and the heUcity correlated ones, both being necessary ingredients of NLO calculations 
employing a subtraction method, have been obtained by properly modifying the routines generated 
by MadGraph. Cross checks have been performed with AutoDipole [28]. Moreover, factorization 
in the soft and coUinear limits has been explicitly checked in double and quadruple precision. 

The finite part of the virtual corrections for /f+jet is extracted from a C++ library of re- 
sults [7, 8], which is compliant with the interface to parton SMC programs proposed in [29]. In 
brief, the virtual amplitudes and loop diagrams associated with eq. (2.1) are all generated analyt- 
ically and, after manipulation with computer algebra programs, automatically translated into C++ 
code. The reduction of tensor integrals up to five-point functions displays good numerical stability. 
Finally, the calculation results in a decomposition of all loop diagrams according to helicity and 
color stracture times scalar functions depending on the extemal momenta only. The infrared (IR)- 
finite scalar integrals are evaluated using the FF package [30, 31] and an efficient caching system 
is appUed to speed-up the entire computation. The current implementation needs approximately 
36 ms CPU time on an Intel Xeon processor with 3.00 GHz for each evaluation of the B function 
(see eq. 4.5 of [13]) at a single phase space point, summing over all contributing subprocesses. 
About half of that time is spent on the computation of the virtual corrections. 

With the input described so far, our POWHEG BOX implementation can be used to generate 
hadronic events accurate to NLO. There is, however, an important technical issue related to multi- 
leg processes which are interfaced to a parton shower: Reactions such as ?f-|-parton possess soft 
and collinear divergences already at LO. Since the POWHEG method relies on generating all events 
starting from a momentum configuration of the underlying Born process, a generation cut, pf'", 
needs to be placed on the transverse momentum of the final-state partons, a procedure well-known 
from standard SMC generators, see e.g., the discussion in [32]. 

This cut for the event generation, which in practice takes values e.g., pf*" Ri 1 GeV, is un- 
physical and has to be always much smaller than the analysis cut on the transverse momentum, 
p™, employed in the definition of the additional jet associated with the tt-pair. The condition 
pT'^ < pT is necessary though not sufficient to provide a suitable sample of unweighted events, 
because the parton shower can in fact increase the transverse momentum in an event generated at 
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Table 1: The dependence of the NLO cross section at Tevatron and LHC on the generation cut p\ and on 
the Bom suppression factor p^"''^ for = = = 174 GeV and CTEQ6M PDFs. Jets are reconstructed 
by the inclusive ^T-algorithm with R=\. 

Bom level with a given < /jf^". This may result in a different number of events that would have 
passed the analysis cut if the generation cut were different. Therefore, it is essential to show that the 
actual dependence on pj™ is negligible for a given fixed analysis cut /?™. We have carefully inves- 
tigated this issue, reporting some examples in Tab. 1 , which shows the independence of the cross 
section on the generation cut pf^'^ for reasonable values of the cut. In particular, this is satisfied for 
/jf*" smaller than pf^/4. Tab. 1 provides also a strong cross check of existing results in [25]. 

Also a second option for dealing with processes with soft and coUinear divergences at Bom 
level is available by generating weighted rather than unweighted events and, thereby suppressing 
the divergences. The POWHEG BOX uses a suppressed cross section B^upp for the generation of the 
underlying Bom configurations [32], 

Bsupp=SxF(pt), (2.3) 

where pj is a measure of the hardness of the extra emission, B denotes the inclusive NLO cross 
section at fixed underlying Bom variables (see eq. 4.2 of [13]), and F{pj) has to be chosen such 
that it ensures finiteness of Bsupp in the limit p^ — 0. Thus, Bsupp becomes integrable and, of course, 
the generated events need to be weighted by l/F(pj). The POWHEG BOX implementation foresees 
the choice 



F{Pt) = 




(2.4) 
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where the variable p-^ controls the extent of the suppression and n is process specific, « = 2 for 
//+jet. In our event generation we have also applied the Bom suppression factor of eq. (2.4) and 
checked its implementation for different p^j^"^ values with very good numerical accuracy, as shown 
in Tab. 1. 

2.1 Checks 

The current implementation has been exposed to a large number of checks. These include, first of 
all, comparisons for a sizable number of differential distributions to NLO with the available fixed 
order results. We have used the settings of [7, 8], i.e., the renormalization and factorization scales 
I^R= I^F =mt = 174 GeV, the analysis cut of p-^ > 20 GeV and p-^ > 50 GeV for the Tevatron and 
LHC configuration, respectively, the parton distribution function (PDF) set CTEQ6M [33], the jet 
algorithm of [34] (inclusive- A:^) with/? = 1 and we have assumed here the top-quarks to be always 
tagged and thus excluded from jet reconstruction. Perfect agreement with [7, 8] has been found 
typically at the per-mille level. Also note that the computation of [7, 8] has employed dipole sub- 
traction with massive partons [35] for the cancellation of the soft and coUinear divergences between 
the real and virtual contributions. As mentioned, the POWHEG BOX relies on FKS subtraction [26] 
and, of course, the results for physical cross sections at NLO must be independent of the chosen 
scheme, which constitutes another strong check of the current implementation. 

The other aspect which demands careful checking is the parton shower. It adds all-order 
perturbative corrections, though to leading logarithmic accuracy only, and the POWHEG approach 
offers one particular choice of how these corrections beyond NLO are included. Since the merging 
with parton showers modifies the QCD predictions, sometimes even for inclusive quantities, one 
should check during the vaUdation process that these modifications reflect real physics effects and 
moreover, are compatible with (un)known higher-order corrections. To that end, we have also 
compared the NLO fixed order and the POWHEG predictions after the first emission, i.e., at the level 
they are written to the Les Houches event file [36] (LHEF from now on), for suitable distributions 
in order to assess similarities and differences between the results. Of particular interest in this 
procedure are those observables, which display great sensitivity to the parton shower effects. Some 
results will be shown below in Sec. 3. 

3. Results 

Throughout this section we present our findings. In comparing to the fixed order NLO results we 
are using the parton shower programs HERWIG [21] (version 6.5.20) and PYTHIA [22,23] (version 
6.4.25), with the default setting of parameters. Moreover, in order to highlight the shower effects 
and to disentangle them from the underlying event (UE) or multiple particle interactions (MPI), we 
have generally switched off these features in the SMC programs. Of course MPI and UE, together 
with the usage of the appropriate SMC tuning, will play an important role when comparing POWHEG 
predictions with data. For the kind of study we are interested in, we have decided not to include 
them in the following plots and we have concentrated on presenting distributions which are not 
extremely sensitive to these effects, where possible. The one exception here is the discussion of the 

'in the POWHEG BOX the value of p^^^ can be assigned at run time trough the bornsuppf act entry of the input card. 
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charge asymmetries, in Sees. 3.1.1 and 3.2.1. There, we will also consider the predictions obtained 
including also MPI and UE in our comparisons in order to asses the stability of our results with 
respect to these effects. 

We have produced different samples of events, to study both, the Tevatron and LHC colUder 
configurations. Each sample contains 20M positive and negative weighted events, produced with- 
out folding (see Appendix A) and with a generation cut pf*° = 2 GeV and 5 GeV, respectively, 
for the Tevatron and the LHC. As already mentioned, the jet reconstruction cut in the analysis has 
been assumed to be > 20 GeV and 50 GeV. We also make use of the bornsuppf act option of 
the POWHEG BOX to damp the low transverse momentum regions in the cross section, artificially in- 
creasing the probabiUty of a hard jet. For this reason, the events have a weight which is the inverse 
of eq. (2.4). The choices p^^^ = 20(100) GeV for the Tevatron (LHC) were adopted. 

For the sake of comparisons and, in order to keep the analysis routine rather insensitive to 
possible contamination in the jet reconstruction procedure, we always force the semi-leptonic decay 
of the (anti-)top-quark when interfacing to a SMC program in the following. Also, as already done 
for the NLO comparisons of Sec. 2, we will always use the jet algorithm of [34] (inclusive-Zcr) with 
R = I and the Sj-recombination scheme.^ In addition, after interfacing to a SMC program and 
performing the hadronization stage, we always exclude from the jet list those jets which happen to 
contain a ^7-flavored hadron, whose origin can be traced back to the (anti-)top-quark decay.^ . 

As we focus entirely on the effect of the parton SMC and the NLO merging, our phenomeno- 
logical analysis in Sees. 3.1 and 3.2 follows [7,8], fixes the renormalization and factorization scales 
Mr =Mf =tnt = 174 GeV, and uses the PDF set CTEQ6M [33]. Studies of the theoretical uncer- 
tainty of the cross section for pp — > tt+ jet due to corrections of higher orders (beyond NLO in 
QCD) have already been conducted in [7, 8] with standard means, i.e., determining the variations 
for the scale choices of /Ir = /ip = mf/2 and = 2mt. In Tab. 2 we report results for the 

Tevatron 1.96 TeV as well as for the 7 TeV and 14 TeV LHC configurations allowing for the inde- 
pendent variations of /7r and /^p of factors m,/2 and 2m, and, at the same time, excluding relative 
ratios of /Vr/^Uf larger than 2 or smaller than 1/2. We find a remarkable small variation, typically 
below 10% . On the basis of the scale variations uncertainty one can thus conclude that higher 
order effects should be small. 

The other main source of theoretical uncertainties, related to the non-perturbative parameters 
such as the PDFs and the associated value of the strong coupling constant as(Mz), have repeatedly 
been addressed in the literature (see e.g., [6]) using modern PDF sets available to NNLO [38-40]. 
For a more recent study concerning the combined uncertainty due to propagated PDF uncertainties 
and uncertainties in a^{Mz) we refer the interested reader to the results discussed in [41]. The top- 
quark mass nit is always taken as a pole mass, and we refer to [42] for a discussion of the running 
mass in the MS -scheme. 

^Comparisons have been done with other jet clustering algorithms, e.g., the anti-fcj- jet algorithm [37] with p-scheme 
recombination, for different values of the parameter R, without significantly changing the general conclusions drawn 
here. 

^This exclusion is applied also when the top-quark decay is performed by POWHEG itself and bottom-quarks coming 
from top-quarks can be present already at the LHEF stage, as it will be the case in Sec 3.3. Other bottom-quarks present 
at partonic level (or B hadron not originating from a (anti-)top-quark) are instead always consistently considered as 
(originating from) light partons and therefore included in the light jets count. 
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pf' [GeV] 


<7nlo [pb] 


TEV 1.96 TeV 
LHC 7 TeV 
LHC 14 TeV 


20 
50 
50 


1.791(2)^^:^?^ 
53.1(2)1^:1 
376.1(2)+20-i 



Table 2: The dependence of the NLO cross section at Tevatron and LHC on the renormaUzation ^ur and 
factorization scales /np. The envelope has been constracted by considering independent variations of /hr and 

/Vp around nit = 174 GeV, by a factor of two in both directions. Combinations resulting in ratios of /Vr//Vf 
larger than 2 or smaller than 1 /2 have been excluded. Jets are reconstructed by the inclusive /ci-algorithm 
with /? = 1 and CTEQ6M PDFs are adopted. 

3.1 Tevatron 

We will begin our studies with the Tevatron (-s/s = 1 .96 TeV), which has provided us already with 
several measurements of differential distributions. These are, in detail, the first measurement of the 
tt cross section dd/dmu, differential in tt-pak invariant mass m^, see [43], the dependence of the tt 
production cross section on the transverse momentum of the top-quark [44] and most prominently, 
the measurements of the forward-backward charge asymmetry in top-quark pair-production [45,46] 
along with its dependence on the ?f-pair invariant mass niti [47,48]. 

In all the plots throughout this article no acceptance cuts are imposed, other than those neces- 
sary to define the hard jet. The lines appearing the legenda are labeled according their origin: NLO 
for the fixed order computation, LHEF for the results after the POWHEG hardest emission without any 
showering and PWG+HER, PWG+PYT for the results obtained by generating the hardest emission with 
the POWHEG method and then performing the remaining shower and hadronization stages with the 
HERWIG or PYTHIA programs, respectively. The lower inset in each plot shows the relative differ- 
ence with respect to the first entry in the legend, which is usually depicted as a solid black line in 
the upper inset. The black dots in the lower inset stand instead as a reference for the zero Une. The 
Monte Carlo integration errors are shown as vertical bars. 

In Fig. 1, we plot the differential cross section as a function of the transverse momentum p^ of 
the top-quark (upper panel) and of the the ff-pair invariant mass mn (lower panel). A first comment 
can be made upon the expected good agreement between the NLO results and the LHEF ones for 
such inclusive quantities. The effect of the shower does not change considerably these results. Only 
at the very low and high ends we observe a significant change which is however still below 30%. 

Similar conclusions can also be drawn looking at the upper panel of Fig. 2, where we plot 
the differential cross sections da/dAi^tt as a function of the azimuthal separation between the top- 
and anti-top-quarks. In particular, this last observable is very stable with respect to the inclusion 
of the parton shower. With exception of the last bin, we observe corrections of a few percent 
only. The understanding is that the parton shower is not able to produce a significant change of 
the relative directions of the two outgoing heavy quarks. In the lower panel instead, we show the 
differential cross sections do /dp" as a function of the transverse momentum of the tt -pair. Shower 
effects at low p" values are now clearly visible. Furthermore, by inspecting the lower inset in the 
plot, showing the relative difference with NLO predictions, it is also possible to see a hint of the 
instabilities that may arise in the fixed order calculation around the pj = 20 GeV jet cut. 
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Figure 1: The differential cross sections as function of the transverse momentum (upper panel) and of 
the f f-pair invariant mass m,f (lower panel), at the Tevatron (y/s = 1 .96 TeV). 
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Figure 2: The differential cross sections as function of the «-pair azunuthal separation A(|),f and transverse 
momentum p", at the Tevatron ~ 1.96 TeV). 
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This behaviour is easily explained considering that the region below the cut may be populated 
only when more than one hard jet is resolved, thus allowing for an imbalance between the hardest jet 
and the tt-pak. Indeed, for the 3-partons ttj configuration, p'^ = by momentum conservation. 
This effectively makes the do/dp!^ distribution for p!^ lower than the jet p-^ cut a LO quantity, 
starting at 0{af). We have verified that the instability gets worse by reducing the bin size and 
that neither the LHEF nor the showered results show a similar behaviour. The same feature will 
also be present in LHC predictions, around the pj = 50 GeV jet cut, as we will show in the lower 
panel of Fig. 7. Such instabilities are well known to arise at any fixed order of the perturbative 
expansion whenever the observable under consideration has a non-smooth behaviour inside the 
physical region or whenever the phase space boundary for a certain number of partons lies inside 
that for a larger number [49]. In our case the LO dc/ p" distribution is discontinuous at the jet p'^"^ 
boundary, since the phase space for 3 partons does not allow p'f = < p^"\ while it is different 
from zero at p!j > It should be stressed that despite these effects are due to unbalanced 

cancellation of higher order soft/coUinear divergences at the critical point(s), they are integrable 
in any finite neighborhood. In ref. [49] it was also demonstrated that the all-order resummation 
of soft divergences restores the infrared finiteness of the predictions also at the critical point by 
the appearance of continuous and infinitely differentiable structures, that were dubbed Sudakov 
shoulders. We have indeed verified that the resummation impUcit in the POWHEG approach does 
reproduce these structures, as can be observed in the figures. 

From Fig. 3 on we investigate the jet structure of the process. To that end, in Fig. 3, we display 
the differential cross sections dcs/dp^"' (upper panel) and dcs/dyjut (lower panel) as function of the 
transverse momentum and rapidity of the hardest jet. Here we notice a moderate shape distortion 
and a different normalization in going from parton level to hadronic events. Both these effects can 
be understood by considering that the jet clustering procedure is significantly affected when going 
from a low multiplicity sample to a fully hadronized event. This is also reflected in the number of 
resolved jets above the minimum pj cut and in their kinematics. 

In the upper panel of Fig. 4 we show the azimuthal difference A(|)jj_^ist between the hardest jet 
and the tt system. A(|)fj_yist is very sensitive to shower effects, since at partonic level the hardest jet is 
constrained to reside in the opposite hemisphere w.r.t. the tt-pair momentum. In the lower panel we 
plot the differential cross section do/dp"^^" as function of the transverse momentum of the system 
composed of the top-quark, the anti-top-quark and the hardest jet. This observable is of particular 
importance, as it is strongly correlated with the transverse momentum of the next-to-hardest jet. 
The latter, in turn, is the quantity where the Sudakov suppression effects should be most evident. 
Due to the jet reconstruction procedure, however, it turns out that those effects are best seen in the 
low p"^'" region. The disagreement in the high Pt^'" tail between NLO and LHEF results, on the 
one side, and showered results, on the other, is instead due to showering and hadronization effects 
and is partly related to our choice of excluding the top-quark from jet reconstruction. A similar 
effect has also been observed (and explained) in single-top production (see Sec. 4 of [50]) and we 
have verified that the same explanation can readily be applied here too, for the tt jut system: At 
high px^'"' there exists an artificial imbalance due to many hadrons coming from the hardest parton 
which are not clustered together with the hardest jet. This creates an effective pj for the system 
containing the hardest jet, giving rise to the harder p"'"" tail observed. 
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Figure 3: The differential cross sections as function of the hardest jet transverse momentum and rapidity at 
the Tevatron = 1.96 TeV). 
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hardest jet (lower panel) at the Tevatron (y/s = 1.96 TeV). 
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Figure 5: The differential cross sections as function of the next-to-hardest jet transverse momentum and 
rapidity at the Tevatron (y^ — 1 .96 TeV). 
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Finally, in Fig. 5 we plot the differential cross section as function of the transverse momentum 
Pt^""* and rapidity yjiai of the next-to-hardest jet. No high-pj enhancement, or suppression, is 
observed in this case. However, large differences in the normalization of distributions are present 
when going from NLO results to LHEF and showered ones. This can again be explained with the 
different number of configurations with are resolved by the reconstruction algorithm as containing 
two hard jets. 

3.1.1 Top-quark forward-backward charge asymmetry 

In view of the measurements of the forward-backward asymmetry in top-quark pair-production at 
the Tevatron [45^7] we have also evaluated the values of the tt rapidity asymmetries in ?f-|-jet 
samples, at various stages of our simulations and with various acceptance cuts. The definitions for 
and A|g used throughout this paper are 




(3.1) 



where Ay,; =yt— yi- Due to the fact that the initial pp state at the Tevatron is a CP eigenstate and 
due to the absence of CP violating effects in QCD, one has ^{yt) = ^X~yt) and Ajfg corresponds 
to Apg in the tt rest frame. When moving to the laboratory frame, where a greater accuracy can be 
reached on the experimental determination of these quantities, Apg gets smaller in magnitude than 
A^g, which in turn is a boost-invariant quantity, depending only on rapidity differences. We report 
both in Tab. 3, expressed as percentage values to ease their readabiUty, at various stages of our 
simulation. We have also included the predictions obtained with the HERWIG and PYTHIA showers 
supplemented with underlying-event and multi-parton interactions activities, respectively, referring 
to them as PWG+HER+UE and PWG+PYT+MPI. In case of PWG+PYT+MPI the Perugia Tune [51] - 
MSTP (5) =320 - has also been adopted. 

For both asymmetries and shower models, we observe that the inclusion of the parton shower 
leads to a significant correction compared to the fixed order and LHEF predictions. For example, 
the results for total asymmetries are almost reduced by a factor of two by the parton-shower. A 
similar effect is also observed in presence of selection cuts on the rapidity or invariant mass of the 
heavy quark pair. This outcome extends the conclusions obtained by previous combined NLO and 
parton shower simulations [47] for the inclusive tt sample."* 

Inspecting the results shown in Tab. 3 more closely, we find that performing PYTHIA and 
HERWIG showers leads to similar predictions for A^: the values are fairly consistent once the 
statistical uncertainties are taken into account. The agreement is instead much worse for Apg. 
The inclusion of underlying-event or multi-particle interaction does not dramatically change this 
picture, the interleaved shower evolution due to MPI and Perugia Tune in PYTHIA results in a 
larger difference with respect to the default behaviour than the simpler UE model in HERWIG. 

The dependence of Afb on the transverse momentum of the «-pair, , has also recently been 
reported in [48]. In ref. [53], the authors have investigated this behaviour in presence of selection 
cuts on the tt pair invariant mass, including also EW corrections. 

^See also [52] for higher-order QCD corrections to ApB and [10] for studies of the tt + 2jet process in particular. 
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Tevatron 1.96 TeV 



NLO [%] 



LHEF [%] 



PWG+HER[%] PWG+HER+UE [%] PWG+PYT [%] PWG+PYT+MPI [%] 



Apg total 



-2.98 ±0.04 -2.95 ±0.05 



-1.75±0.11 



-1.70±0.11 



-1.49±0.11 



-1.36±0.11 



A'pji, \yt\<^-0 
A^B> bfl>i-o 



-2.60 ±0.04 
-6.38±0.19 



-2.55 ±0.05 
-6.51 ±0.15 



-1.51±0.12 
-3.79 ±0.35 



-1.53±0.11 
-3.15±0.34 



-1.31±0.12 
-2.99 ±0.35 



-1.22±0.12 
-2.58 ±0.34 



Al,n , \mtt \ < 450 GeV 



^FB 5 
^FB ' 



AU , \mtt \ > 450 GeV 



-1.90 ±0.06 
-4.70 ±0.06 



-1.80 ±0.06 
-4.77 ±0.08 



-1.24±0.14 
-2.70±0.19 



-1.24±0.13 
-2.54±0.18 



-0.81 ±0.14 
-2.66±0.18 



■1.00±0.14 
-1.98±0.18 



A^B> p¥> lOGeV 
A^B ' p't ^ 20 GeV 
A^B > Pt > 35 GeV 
A^B ' Pt > 50 GeV 
A^B ' Pt > 75 GeV 



-2.95 ±0.04 
-2.41 ±0.05 
-3.90 ±0.06 
-4.31 ±0.07 
-4.88 ±0.08 



-2.93 ±0.05 
-2.94 ±0.05 
-3.85 ±0.05 
-4.33 ±0.06 
-4.62 ±0.08 



-2.64 ±0.06 
-2.80 ±0.05 
-3.54 ±0.06 
-4.00 ±0.07 
-4.33 ±0.09 



-2.59 ±0.06 
-2.80 ±0.05 
-3.55 ±0.06 
-4.02 ±0.07 
-4.29 ±0.09 



-2.58 ±0.06 
-2.85 ±0.05 
-3.67 ±0.06 
-4.19 ±0.07 
-4.59 ±0.09 



-2.39 ±0.06 
-2.81 ±0.05 
-3.63 ±0.06 
-4.19 ±0.07 
-4.56 ±0.09 



AJJb total 


-4.40 ±0.04 


-4.34 ±0.05 


-2.80 ±0.11 


-2.54 ±0.11 


-2.22 ±0.11 


-1.84±0.11 


A«B, \Ayn\<LO 
A^B, |Am|>1.0 


-2.70 ±0.04 
-19.48±0.18 


-2.62 ±0.05 
-19.54 ±0.22 


-1.71±0.11 
-10.52 ±0.52 


-1.91±0.11 
-9.75 ±0.51 


-1.39±0.11 
-9.22 ±0.52 


-1.16±0.11 
-7.54±0.51 


ApB ' I'^wl < 450 GeV 
A^B , > 450 GeV 


-3.59 ±0.06 
-5.70 ±0.06 


-3.51 ±0.06 
-5.66 ±0.08 


-2.67 ±0.14 
-3.03±0.19 


-2.36±0.13 
-2.88±0.18 


-1.74±0.14 
-3.06±0.18 


-1.63±0.14 
-2.20±0.18 


A^B , Pt > 10 GeV 
A^B ' Pt > 20 GeV 
A§s , Pt > 35 GeV 
A^B , Pt > 50 GeV 
Afs , p'^ > 75 GeV 


-4.35 ±0.04 
-3.71 ±0.05 
-5.72 ±0.06 
-6.25 ±0.07 
-6.62 ±0.08 


-4.32 ±0.05 
-4.29 ±0.05 
-5.52 ±0.05 
-6.11 ±0.06 
-6.45 ±0.08 


-3.98 ±0.06 
-4.22 ±0.05 
-5.16±0.06 
-5.70 ±0.07 
-5.99 ±0.09 


-3.86 ±0.06 
-4.18 ±0.05 
-5.17 ±0.06 
-5.74 ±0.07 
-5.94 ±0.09 


-3.72 ±0.06 
-4.15 ±0.05 
-5.21 ±0.06 
-5.85 ±0.07 
-6.27 ±0.09 


-3.51 ±0.06 
-4.11 ±0.05 
-5.21 ±0.06 
-5.92 ±0.07 
-6.25 ±0.09 



Table 3: Results for the forward-backward asymmetries Apg and Apg at the Tevatron, for the tt j sample with various acceptance cuts. Jets are reconstructed by 
the inclusive ^T-algorithm with R= I, above the p^^ > 20 GeV minimum jet cut. 



Since this observable only receives contributions starting at O {(xl), even for the tt inclusive 
sample - an extra jet is always required to have a non-zero p'.J - our approach offers the possibihty 
to evaluate it at NLO supplemented with shower effects for the first time. This should improve over 
the MCSNLO and POWHEG simulations foTtt production, that are only leading-order/leading-logarithm 
accurate for Afb when p'^ > 0. 

In Tab. 3 we have also included our predictions for the forward-backward asymmetries with 
p'^ > 0. We now see, that the results after showering with two different SMCs are consistent almost 
everywhere, apart from the very low p'J region, which thus contributes to generate the differences 
in A^g. Moreover, for both asymmetries Apg and A^g, we observe that for Pj> 10 GeV the 
inclusion of the parton-shower leads to results in fair agreement with the NLO ones. The large 
differences between NLO and showered results observed in the total asymmetries are thus traced 
back to the region of small p'^'s. This appears plausible since in that region soft and colUnear 
emission produced by the parton shower will be more important compared to the large p'-l region. 
We also note that our findings are different from the observations made in [48] where a stronger 
effect of the parton shower using specific tunes was observed over the entire range. One has 
to remind, however, that for a leading-order SMC such as the PYTHIA program used in [48], the 
contributions to the asymmetry for pfj > are generated exclusively by the shower. It thus looks 
feasible that a change in the parton shower will be reflected in a overall change of the predictions. 

Whether the aforementioned large corrections produced by the parton shower could be re- 
sponsible in part for the observed discrepancy between theory and experiment is however unclear 
because in the tt sample this effect may cancel. Given the stability of the theoretical prediction 
above p'-} = 10 GeV — now available at NLO accuracy — a detailed comparison of experimen- 
tal data and theory predictions may provide useful information to further scrutinize the observed 
discrepancy in the inclusive sample. 



3.2 LHC 

In Figs. 6 to 11 we report a similar set of distributions as already shown for the Tevatron, this 
time for the LHC colUder at ^/s = 1 TeV. For most of the observables, the same comments and 
conclusions as reported for the Tevatron case are, in general, vaUd here, too. For this reason, we 
will not repeat the same remarks here. There are, however, some new plots that deserve discussion 
and explanations. 

In the lower panel of Fig. 6 we show the differential cross section as a function of the rapidity 
distribution of the top-quark jt. This inclusive quantity shows the expected good agreement be- 
tween NLO and LHEF results. It is also pretty stable with respect to the inclusion of parton shower 
effects, resulting only in a few percent change in the overall normaUzation. 

In the upper panel of Fig. 9 we show the differential cross section do/dp^ '^i as function of the 
scalar sum of the relative transverse momenta of the particles in the first jet, /Jx'reb which is defined 
with respect to the jet axis in the frame where the jet has zero rapidity. In general, we define for the 
n-th jet, 

Mb _ y 1^' ^^^""1 (2 2) 
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where ki denotes the momentum of the / particle in the n-th jet, see also [50]. The observable 
behaves as expected: The NLO curve shows the unphysical enhancement associated with the IR 
singularity at p{^^l^ = 0, the LHEF result corrects for this only partially, owing to the Sudakov sup- 
pression by final-state radiation, while the fully showered events manifest the physical suppression 
at low pxlrei. Moreover, in the lower panel of Fig. 9 we show the distribution of the invariant mass 
of the system made by the top-quark and the hardest light jet. 

Eventually, in the lower panel of Fig. 10, the total exclusive cross section as a function of the 
number of resolved jets, Njets, is shown. The moderate difference in the Njets = 2 bin between 
the NLO and LHEF results may be considered a consequence of the different normalization of the 
distributions involving the next-to-hardest jet, as shown e.g., in Fig. 11. Higher bins are instead 
populated only by the shower and thus predictions for those bins are only accurate in the strict 
soft/collinear limit. 

One noticeable difference w.r.t. Tevatron results is the reduced discrepancy between the show- 
ered and partonic results for observables that involve the hardest or next-to-hardest jets, e.g.. Figs. 8 
and 1 1 . A possible explanation for this is in the higher jet cut (px = 50 GeV) used for LHC analysis, 
that largely restricts the regions where shower effects may become important. 

A first phenomenological application of the POWHEG implementation discussed here has already 
been reported in [54], namely the study of ttj-events at the LHC as a function of the invariant mass 
of the multi-jet systems as a means of measuring the top-quark mass. 

3.2.1 Top-quark charge asymmetry 

In Tab. 4 we present our estimate of the charge asymmetries in the tt j sample, measured with 
respect to the (anti-) top-quark pseudo-rapidity r\ and rapidity y, at various stages of our simulation. 
The two observables are defined in the same way 



with A|x| = 1^:, I — \xf \ and x = r\ or x = y, respectively. By looking at our results, one can see 
that a positive asymmetry is present for both observables, albeit larger in magnitude for A^. More- 
over, results obtained with the inclusion of PYTHIA and HERWIG showers are in fair agreement, even 
after the addition of MPI or UE. We also observe that the numerical evaluation of the asymmetry is 
rather challenging at fixed order, leading to large statistical uncertainties. Only after increasing the 
statistics by one order of magnitude we were able to obtain errors comparable with those obtained 
after the inclusion of the parton shower. We consider the results including the effects of the parton 
shower as the best prediction currently available for the asymmetry. We refer to refs. [55, 56] for 
recent LHC measurements of A^ and A^. in the tt inclusive sample. 

In closing we remark, that we have also carefully tested the LHC nominal center-of-mass 
energy, producing the same set of distributions as in Figs. 6 to 1 1 for -^/j = 14 TeV collisions. 
Since no significant difference in the behaviour has been observed in any of the distributions, with 
respect to the 7 TeV case, we refrain from including the corresponding plots here. 
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Figure 6: The differential cross sections as function of the transverse momentum and of the rapidity yt 
of the top-quark at the LHC (^/s = ^ TeV). 
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Figure 7: The differential cross sections as function of the f f-pair invariant mass niif and transverse momen- 
tum pt} at the LHC {^/s = l TeV). 
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Figure 8: The differential cross sections as function of the hardest jet transverse momentum and rapidity at 
the LHC (V? = 7 TeV). 
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Figure 9: The differential cross sections as function of the scalar sum of transverse momenta inside the 
hardest jet (upper panel) , as defined in eq. (3.2), and the invariant mass of the system made by the top-quark 
and the hardest jet (lower panel), at the LHC (y/s = 7 TeV). 
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LHC (^A = 7 TeV). 



-22- 



> 

0) 

O 

\ 



b 
< 



ft 



\ 
b 



b 
< 



10 







— , 10 



-1 _ 



10 



-2 _ 



i!; 10-3 



b 

lo-'^ 



10 



-5 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




NL^ 

+- LHEF 
■*--PWG+HER 
■®--PWG+PYT 



LHC 7 TeV 
ki— algo R=l 
p^** > 50 GeV 

I .... I 




J I I I I I I I I I I I I I I I I I I I I L 



0.5 
0.0 
-0.5 



T 



i. . I 



I ' ' ' ' I 



J I I I I I I I I I I I I I I I I i_ 



100 200 300 



400 



500 



10.0 
5.0 

2.0 
1.0 
0.5 

0.2 I— 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



LO 
LHEF 
-■»*--PWG+HER — I 
PWG+PYT 



0.1 

0.4 

0.2 

0.0 

-0.2 
-4 



LHC 7 TeV 
kT-algo R=l 
> 50 GeV 



J I I I I I I I I I I I I I I I I I L 



1 



zn — I — -T — I — |- 



-1 — I — I — r 



1^ 1 1 1 1 — -j 1 1 1 ITZ 

J I I I I I I I I I I I I I I I I I I : 



-2 



Figure 11: The differential cross sections as function of the next-to-hardest jet transverse momentum and 
rapidity at the LHC (V^ = 7 TeV). 
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LHCTTeV NLO [%] LHEF [%] PWG+HER [%] PWG+HER+UE [%] PWG+PYT[%] PWG+PYT+MPI [%] 



0.19 ±0.09 0.18 ±0.06 0.46 ±0.10 0.26 ±0.11 0.40 ±0.11 0.57 ±0.11 

Al 0.51 ±0.09 0.47±0.06 0.73±0.10 0.52±0.11 0.66±0.11 0.76 + 0.11 



Table 4: Results for the top-quark charge asymmetry Ac in the tt j sample at the LHC (-y/i = 7 TeV). Jets are reconstructed by the inclusive ^T-algorithm with 
R=\, above the pj^ > 50 GeV minimum jet cut. 



Tevatron 1.96 TeV LHEF [%] PWG+HER [%] PWG+HER+UE [%] PWG+PYT [%] PWG+PYT+MPI [%] 



A*,B -1.83±0.05 -0.94±0.11 -0.92±0.11 -0.57±0.11 -0.41 ±0.11 

A^^~ -2.21 ±0.05 -1.00 ±0.11 -0.99 ±0.11 -0.59 ±0.11 -0.46 ±0.11 



PWG+HER' [%] PWG+HER'+UE [%] PWG+PYT' [%] PWG+PYT'+MPI[%] 



0.69±0.11 -0.93±0.11 -0.71 ±0.11 -0.47±0.11 

1.13±0.11 -1.11±0.11 -0.90±0.11 -0.85±0.11 



Table 5: Results for the lepton charge asymmetries Apg and A^g in the tt j sample at the Tevatron {^/s = 1.96 TeV). Jets are reconstructed by the inclusive 
^T-algorithm with R=l, above the p^^ > 20 GeV minimum jet cut. 



A^ 

^FB 

^FB 



3.3 Top-quark decay 

The studies so far have been based on the assumption of stable (anti-)top-quarks during both the 
fixed order calculation stage and the generation of the hardest radiation according to the POWHEG 
method. Thus, the output resulting from these stages includes on-shell top-quarks and extra light 
partons only. When interfacing to a SMC program for performing the rest of the showering and the 
hadronization steps, no information concerning the spin or helicity of these top-quarks is retained 
and the decay is performed by the SMC, which averages over the possible top polarizations. 

In this way one does not correctly account for top-quark off-shell effects, non-resonant produc- 
tion mechanisms and for the top-quark spin correlations. However, given the large ratio between 
the top-quark mass and the top-quark width, the first two effects are suppressed as 0{rt/mt) and, 
therefore, are rather modest. This was also explicitly shown for the tt production case by the recent 
calculation of the NLO corrections to the W~^W~bb final state [57,58]. 

Top-quark spin correlations at hadron colliders were instead realized to be an important tool 
to study top-quark properties since the early work of [59]. Also, the QCD corrected spin analyzing 
power of jets in decays of polarized top-quarks has been extensively discussed in [60]. The possi- 
bility to measure these correlation effects is particular to top-quarks because of their large mass and 
short lifetime, which prevent non-perturbative QCD effects from significantly depolarizing them 
before the decay. 

The usual approach to study this phenomenon relies on taking the top-quark zero-width limit 
and separately considering NLO corrections to the production and decay processes. This is done by 
means of a spin correlation matrix which accounts for all the possible spin configurations occurring 
in the two stages. Several studies following this approach have akeady been presented, concerning 
both tt production [59] with spin correlations at NLO and the tt j case [20] with spin correlations 
accounted for at leading order. 

Here we follow a different approach, that was first introduced in [61] and has already been 
employed and discussed in both implementations, POWHEG [18, 19, 62] and MC@NLO [16, 63]. In this 
way we avoid having to introduce the spin correlation matrices, but we nevertheless recover the 
exact leading-order spin correlations. To be more precise, we have LO spin-correlations in the soft 
and coIUnear regions, while the same accuracy of the complete NLO calculation, namely the matrix 
element, 

PP{PP) ^{t^ M+v) {t ^ brv)jj , (3.4) 

is used for the hard regions.^ For the sake of brevity, we explain this method here only briefly, 
and instead, refer to the aforementioned publications for more details. The basic idea is to first 
generate events with stable top-quarks (un-decayed events) through the usual POWHEG machinery 
and then generate the decay products according to the matrix element for the full production and 
decay process (decayed events). 

Moreover, in order to produce more realistic final states, we have also allowed for a reshuffling 
of the momenta of the top-quark decay products, resulting in off-shell top-quarks and W bosons, 
whose virtuaUties have been distributed according to Breit-Wigner shapes. This reshuffling is such 

^In case no radiation harder than = 0.8 GeV is generated by the POWHEG method, the event is classified as 
Born-like: the corresponding matrix elements for the decayed process , i.e., pp{pp) {t M+v)(f M~v)j, are 
then employed. 
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that the relative 3-momentum of the (anti-)top-quark in the tt rest frame and the 3-momentum of the 
radiated Ught partons in the partonic center-of-mass system are kept fixed. It also takes care of the 
small changes in the phase space and luminosity due to the t and t quark being slightly off-shell with 
respect to the chosen pole mass value. In any case, we restrict the top-quark and W virtualities off- 
shellness to be in between 10 widths. We have checked that this modification does not significantly 
change the results obtained in the top-quark and W boson zero- width approximation, by comparing 
several distributions for the two cases.^ 

Furthermore, since charged leptons (or down-type quarks) coming from the top-quark decays 
are the most sensitive probe of the top-quark spin directions, we concentrate on the di-leptons 
chaimel only. We thus always assume t W~^b — > i~^vb and t W~b — > i~Vib, with branching 
ratio BR{W £v) = 0.108. However, we stress that our implementation is completely general for 
what concerns the top-quark decay products, such that the lepton-t-jets or the all hadronic search 
channels can be simulated as well, if desired.' 

Thus, taking as an example the semi-leptonic top-quark decay channel, the conditional proba- 
bility for the decayed events, starting from un-decayed ones, is obtained from 

dP(^. 10) , ) = i ^ec. (<^>undec. , ^,^Mv,^J^btd .^^^ 

BR{t^bev)BR{t^bfv) ^Wundec.(*J>u„clec.) 

where Oj^^. = {<I>undcc 5*I'r^Mv,*J*r^Mv} is the full space phase including the decay products, while !M's 
are the matrix elements squared for the given partonic process. In practice, however, to generate 
^>dec. configurations according to eq. (3.5), one makes use of an upper bound ?7dec. for the ratio 
^ec./^^^^indec. and then proceed as follows: 

I. Generate a tentative decay kinematics Ojec, starting from <I>u„dK.- 

n. Extract a random value r in the range [0, f/dec.]- 

in. If r < ^Hiec.(*J*dec )/^Wundec.(*J*undec.)> then the decay kinematics is allowed and the event is 
accepted. Otherwise go back to step I. 

An efficient functional form of the upper-bounding function may be readily reconstructed from 
the structure of top-quarks decay: 



t^dec. ^dec. = ^ 7-2 1^1 7771 / 2 ^2 X ^ ^ > (3-6) 

{Mf-mfY + miTf {M^^- m^f + m^T^ 



where Mj and = [ki+ky)^ are the top-quark and W boson virtualities and fMt^bw and M^^j^ 
are the squared amplitudes for the 1 — > 2 decay processes. The choice of the normalization factor 
fA^ has been performed by sampling the allowed phase space, by requiring that the inequality 

iHl{c.(*^*u„dec.,<J>dec.) < ^K^dec. (*^*undec. ) t^dec. (<I>dec. ) (3.7) 

always holds, for any given specific subprocess flavour /. 



^The zero-width approximation in the decay may always be enforced by setting the zerowidth token to 1 in the 
input card. 

'See e.g., [64] for a review of observables sensitive to spin-correlation effects in the lepton+jets channel. 
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In the rest of this section we present results obtained within this approach. We have set 



n = 1.31 GeV , mw = 80.398 GeV , Tw = 2.141 GeV , 

ttEM = 1/128.89 , sin^ 6^ = 0.22265 , (3.8) 

and we have considered a simplified form for the CKM matrix, with a mixing between the first two 
generations only: 

K^ = K. = 0.975, Vus = Vcd = 0.222, V.^, = 1 (3.9) 

and all the other entries set to zero. No extra acceptance cuts on the leptons are placed in the 
following plots. 

With these ingredients we are able to offer predictions for top-quarks decay products, at the 
LHEF level and after shower and hadronization. We compare them in Fig. 12 , where we plot the 
positively charged lepton transverse momentum and the rapidity of the negatively charged one, at 
the 7 TeV LHC and 1 .96 TeV Tevatron colUders, respectively. We note that in these distributions 
shower effects are rather small. Moreover, one does not expect to see any dependence on the 
top-quark spin. 

Spin correlation effects in tt production are usually studied by choosing t and t spin quantiza- 
tion axis and looking at reference frames and angular distributions that are particularly sensitive to 
these correlations^, a procedure that can easily be performed in SMC simulations. However, due 
to the presence of unobserved neutrinos, this can become a difficult task in the data analysis for an 
experimental collaboration. An example of such observables is the double differential distribution, 

1 d^a 

a (icos0i<icos02 ' 

where 0] and 02 are the angles between the directions of the flights of the leptons coming from the 
decayed top-quark in the t (t) rest frame and an arbitrary direction that defines the quantization axis 
for the (anti-)top-quark spin. 

The double differential distribution in eq. (3.10) can be evaluated to good approximation to 

^(l-KcosBicosei) , (3.11) 

which is only vaUd, though, for fully inclusive cross sections. As soon as cuts are imposed the form 
of the distribution will change and a further dependence on cos 612 could be introduced.^ The cuts 

may also affect the simple interpretation of K, which in general depends also on the choice of the 
spin quantization axis. Nevertheless even in the presence of cuts the distribution is still a useful 
observable to study the impact of spin correlations. This quantity has been recently measured by 
DO [66] for the choice of the beam direction as quantization axis, finding K = 0.10 it 0.45. Due 
to the large uncertainties, this value is both compatible with no spin correlation and with the NLO 
predictions K = 0.777 of [59]. 



^See e.g., [65] for strategies to maximize the spin correlation of ff -pairs produced at tiie LHC. 

'We note that electroweak effects and QCD absorptive parts lead to a tiny polarization of the top-quark and thus 
a deviation from the simple form in eq. (3.11) of the double differential distribution shown in eq. (3.10), even in the 
absence of cuts. 
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Figure 12: The differential cross sections as function of the positively charged lepton transverse momenta 
at the LHC (y/s = 7 TeV) and the the differential cross sections as function of the negatively charged lepton 
rapidity at the Tevatron (y/s = 1 .96 TeV). 



-28- 



0.100 - 



-g, 0.050 



CM 

q> 

GO 

o 
o 



m 
O 
o 



b 

X) 



b 
< 



OS 

m 
O 
o 



b 
T3 



b 
< 



0.020 h 
0.010 
0.005 

0.002 I— 



-1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



— PWG+HER unc. 
- PWG+HER corr 



_i I I I I I I i_ 




TEV 1.96 TeV 
kj— algo R=l 
> 20 GeV 



_1 I I L 



0.5 = 

0.0 
-0.5 H 



-1 — I — I — I — I — I — I — I — r 



_i I I I I I I I I I I I I I I I I I i_ 



-1.0 
0.100 



0.050 



CM 
CD 

w 0.020 

o 

o 



-0.5 0.0 

cos 6i cos 02 



0.5 



1.0 



0.010 
0.005 

0.002 I— 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



PWG+PYT unc. 

- PWG+PYT corr.— I 




TEV 1.96 TeV 
k^— algo R=l 
p^'* > 20 GeV 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

0.51 I I I 

0.0 
-0.5 



1 



1 



1 



-1.0 -0.5 0.0 0.5 

cos di COS 02 



1.0 



Figure 13: Effect of the inclusion of spin correlations when interfacing to HERWIG (upper) and PYTHIA 
(lower) programs. Results for the Tevatron collider (^/s = 1.96 TeV). 
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Figure 14: Effect of the inclusion of spin correlations when interfacing to HERWIG (upper) and PYTHIA 
(lower) programs. Results for the LHC collider (yG ^7 TeV). 
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In Fig. 13 we present our predictions for the distribution in eq. (3.10) at the Tevatron using 
the beam-axis to define the lepton angles, obtained by including the spin correlation effects as 
explained above in our tt j sample. We compare these results with those obtained by letting the 
SMC program performing the uncorrected top-quark decays. The effect of the spin correlations is 
clearly visible, leading to differences as large as 50% at extreme cosGi cos 62 values. 

Alternatively, one can look at the distributions of the top-quark decay products in the labora- 
tory frame. In this case the problem is to devise clever enough cuts to single out spin-correlation 
effects which are usually hardly visible.'" 

There are, however, notable exceptions. One such example, the azimuthal distance between 
the leptons coming from the top- and the anti-top-quark, has already been shown in [70] to depend 
mildly on spin-correlations. Results for this observable has also been presented in [71], for the on- 
going LHC 7 TeV run. Later on, it has been suggested that imposing the mft < 400 GeV acceptance 
cut should enhance this dependence in the di-lepton search channel at the LHC [64]. However, due 
to the inherent ambiguities of the unobserved neutrinos, the true md value cannot be reconstructed 
on an event-by-event basis in a real experiment. In our case, we simply looked back in the shower 
history to reconstruct unambiguously the true m^. A more realistic study would use the definition 
of ruti value averaged over all the possible neutrino assignments proposed by the authors of [64]. 
However, as shown in [64], this only leads to minor differences. 

In Fig. 14 we show our results with both the HERWIG and P YTHIA showers and the hadronization 
stage simulations. In both panels, spin correlation effects are again clearly visible as a change in 
the slope of the distribution. 

As a final remark, we investigate the effects induced by top-quark charge asymmetries on 
the di-lepton final state, by studying the leptonic charge asymmetries. This issue has also been 
addressed in a recent paper [72], for the « inclusive production, including next-to-leading order 
QCD and EW corrections. In ref. [20], the authors provide predictions for this quantity at NLO in 
the tf j sample. The definitions for the asymmetries follow the ones in eq. (3.1), but are expressed 
in term of leptonic rapidities in the laboratory frame: 




(3.12) 



with Ayi+£- = yi+ —yi-. We present our results in Tab. 5, also reporting, for ease of compar- 
isons, the results including UE and MPI and the ones for the uncorrelated case, which we dub 
PWG+HER' and PWG-i-PYT', where the (anti-)top-quark decay is performed by the SMC neglect- 
ing spin-correlation effects . Comparing the results including the spin correlation to the results 
where the top-quark decay is handled by the shower program — spin correlations are thus not 
taken into account in the latter predictions — we observe only a small difference between the two 
approaches. On the other hand a significant change can be observed when comparing LHEF re- 
sults with the results including shower effects. Our conclusion is that most of the effect present at 
parton level is washed out by the parton shower. Indeed, investigating the behaviour of the leptonic 

1" A recent paper [67] suggests the use of these distributions to construct a likelihood function for the spin-correlation 
hypothesis, that can be then tested statistically for its validity. See also [68,69] for up-to-date results using this approach 

at Tevatron. 
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forward-backward charge asymmetries in presence of a cut on the transverse momentum of the in- 
termediate tt pair, one can trace back the large differences between results at partonic and hadronic 
level to the region of low pif 's once again, as it was the case for the top-quark forward-backward 
charge asymmetry shown in Tab. 3. 

4. Conclusions 

We have studied top-quark hadroproduction in association with a jet merged with parton showers 
to NLO in QCD. For that purpose, we have presented an implementation of the process ff+jet in 
the framework of POWHEG BOX. A careful validation has been performed by comparing with known 
results. For Tevatron and LHC at = 7 TeV we have investigated a large number of distributions, 
some of them with particular sensitivity to parton shower effects, such as the differential cross sec- 
tions as function of the transverse momentum of the ?f-pair and the hardest jet. We have presented 
a detailed study of the impact of the parton-shower on the top-quark charge asymmetry. We find 
that the inclusion of the shower changes significantly the NLO predictions. A detailed analysis 
allowed us to trace back the origin of these corrections to the low pj region. We stress that the 
results presented here represent the first NLO evaluation of this distribution. Excluding the region 
p'l <\Q GeV we find that the parton shower leads only to a marginal change of the charge asym- 
metry biimed \a p^^. A detailed comparison with future measurements may help in clarifying the 
discrepancy observed in the inclusive forward-backward charge asymmetry. However, we stress 
that in this work, we always refer to the asymmetries defined for the tt j sub-sample of the inclusive 
tt production. 

We have also taken into account top-quark decays at LO in QCD, so that spin correlations, 
e.g., for the leptons from the decaying «-pair can be studied. Similar to the case of inclusive 
top-quark pair production the effects of spin correlations are clearly visible in double differential 
distributions of the top-quark (and anti-quark) decay products. Since tt+']Q\, events represent an 
important fraction of all tt events, the measurement of this effect should be feasible at Tevatron and 
LHC. 

The present work offers several options for extension. For instance, we have restricted our- 
selves in this article entirely to the study of top-quark hadroproduction, although the available QCD 
corrections at NLO [7, 8] are equally appUcable to the process ^7^-|-jet with massive bottom-quarks. 
Therefore, the investigation of ^7^-|-jet in the framework of POWHEG BOX is feasible, even if further 
refinements in the code will be required to perform this task. 

Most important, however, will be an extension of the phenomenological analysis presented in 
this paper to meet the real conditions of the experimental environment, e.g. , studying more carefully 
tuning, UE and MPl effects of the SMC and, potentially, for other acceptance cuts. Also, results 
with different values of the top-quark mass, strong coupling constant as(Mz), n^, I^y scales and 
with modem sets of PDFs will be needed. 

Another natural development is the application of the merging procedure recently proposed 
in [73] to the MENLOPS improved tt and ttj samples, in order to improve the description of both 
fully inclusive // and inclusive tt+]ei productions. However, this will be the subject of a separate 
publication. 
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In general, we believe that a more thoroughly phenomenological study would be better per- 
formed in the framework of an experimental collaboration. Therefore, in order to facilitate all 
these future studies, the code of the POWHEG BOX implementation for the process tt+jet will be 
made publicly available on the web-page http : //powhegbox .mib . inf n . it. 
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A. Folding 

In this appendix we briefly review the folding procedure, its implementation in the POWHEG BOX 
and the application to the present calculation. For a more detailed explanation of the procedure 
itself we refer to [32]. 



folding : 




1-1-1 


5-1-1 


1-5-1 


1-1-5 


5-5-5 


5-5-10 


TeV 1.96TeV,;?f" 


= 2GeV 


0.20 


0.17 


0.14 


0.16 


0.07 


0.06 


LHC 7 TeV, pf'^" = 


= 5GeV 


0.22 


0.19 


0.16 


0.19 


0.12 


0.11 


LHC 14 TeV, /rf" : 


= 5 GeV 


0.28 


0.23 


0.26 


0.19 


0.17 


0.16 



Table 6: The fraction negative weighted events generated using different folding factors. 

The POWHEG method allows for the generation of positive weighted events only. This is a 
direct consequence of the perturbative nature of the exclusive differential cross section integrated 
over the Bom variables B (see e.g., eq. 4.2 in Sec. 4 of [13] and comments thereafter). In our 
implementations this integration is performed in the Monte Carlo approach: we generate fully 
exclusive configurations which depend also on the three variables associated to the extra radiation 
(^,y,(^ in the FKS approach [26]) and then we sum over them discarding the extra dependency. 
However, the weight of this fully exclusive configurations, which is dubbed B, is not guaranteed 
to be always positive, in its enlarged phase space. The positivity of B may always be recovered 
by keeping the Bom variables fixed and redefining B as the average over more phase space points, 
which differs only in the radiative variables. This clearly corresponds again to integrating over the 
radiative variables. This is the essence of the folding procedure: the number of phase space points 
considered in each evaluation of the B, for each radiative variable, is called folding factor. 

An immediate drawback is that the folding procedure results in a increase of the program run 
time which is proportional to the product of the folding factors for the three radiation variables. In 
particular cases, it tums out to be more practical including negative weighted events in the final 
sample. However, these performance costs in the generation stage may be well balanced if full 
detector simulations are included or if it is required to have positive weights only. For this reasons, 
we report in Tab. 6 the observed percent occurrence of negative weighted events in case of different 
folding levels and different colliders. 



-33- 



References 

[1] P. Nason, S. Dawson, and R. K. Ellis, The one particle inclusive differential cross-section for heavy 
quark production in hadronic collisions, Nucl. Phys. B327 (1989) 49-92. 

[2] W. Beenakker, W. L. van Neerven, R. Meng, G. A. Schuler, and J. Smith, QCD corrections to heavy 
quark production in hadron-hadron collisions, Nucl. Phys. B351 (1991) 507-560. 

[3] M. L. Mangano, R Nason, and G. Ridolfi, Heavy quark correlations in hadron collisions at next-to- 
leading order, Nucl. Phys. B373 (1992) 295-345. 

[4] M. Czakon and A. Mitov, Inclusive heavy-flavor hadroproduction in NLO QCD: the exact analytic 
result, Nucl.Phys. B824 (2010) 111-135, [0811 .4119]. 

[5] R. Bonciani, S. Catani, M. L. Mangano, and R Nason, NLL resummation of the heavy-quark 
hadroproduction cross- section, Nucl. Phys. B529 (1998) 424-450, [hep-ph/9801375]. 

[6] M. Aliev, H. Lacker, U. Langenfeld, S. Moch, P. Uwer, et ah, HATHOR: HAdronic Top and Heavy 
quarks crOss section calculatoR, Comput.Phys.Commun. 182 (2011) 1034-1046, [1007 . 1327]. 

[7] S. Dittmaier, P. Uwer, and S. Weinzierl, NLO QCD corrections to tt + jet production at hadron 
colliders, Phys.Rev.Lett. 98 (2007) 262002, [hep-ph/0703120]. 

[8] S. Dittmaier, P. Uwer, and S. Weinzierl, Hadronic top-quark pair production in association with a 
hard jet at next-to-leading order QCD: Phenomenological studies for the Tevatron and the LHC, 
Eur.Phys.J. C59 (2009) 625-646, [0810 . 0452]. 

[9] G. Bevilacqua, M. Czakon, C. Papadopoulos, and M. Worek, Dominant QCD backgrounds in Higgs 
boson analyses at the LHC: a study of pp ^ tt -\-2 jets at next-to-leading order, Phys.Rev.Lett. 104 
(2010) 162002, [1002.4009]. 

[10] G. Bevilacqua, M. Czakon, C. Papadopoulos, and M. Worek, Hadronic top-quark pair production in 
association with two jets at Next-to-Leading Order QCD, 1108.2851. 

[11] S. Frixione and B. R. Webber, Matching NLO QCD computations and parton shower simulations, 
JHEP 06 (2002) 029, [hep-ph/ 2 4 2 4 4 ] . 

[12] P. Nason, A new method for combining NLO QCD with shower Monte Carlo algorithms, JHEP 11 
(2004) 040, [hep-ph/0409146]. 

[13] S. Alioli, P. Nason, C. Oleari, and E. Re, A general framework for implementing NLO calculations in 
shower Monte Carlo programs: the POWHEG BOX, JHEP 06 (2010) 043, [1002 .2581]. 

[14] R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, R. Pittau, et al. Scalar and pseudoscalar Higgs 
production in association with a tt-pair, Phys. Lett. B701 (2011) 427^33, [1104 .5613]. 

[15] S. Frixione, P. Nason, and B. R. Webber, Matching NLO QCD and parton showers in heavy flavour 
production, JHEP 08 (2003) 007, [hep-ph/0305252]. 

[16] S. Frixione, E. Laenen, P. Motylinski, and B. R. Webber, Single-top production in MC@NLO, JHEP 
0603 (2006) 092, [hep-ph/0512250]. 

[17] S. Frixione, P. Nason, and G. Ridolfi, A Positive-Weight Next-to-Leading-Order Monte Carlo for 
Heavy Flavour Hadroproduction, JHEP 09 (2007) 126, [0707 . 3088]. 

[18] S. Alioli, P. Nason, C. Oleari, and E. Re, NLO single-top production matched with shower in 
POWHEG: s- and t-channel contributions, JHEP 09 (2009) 1 1 1 , [0 9 7 . 4 7 6] . 



-34- 



[19] E. Re, Single-top Wt-channel production matched with parton showers using the POWHEG method, 
Eur.Phys.J. C71 (2011) 1547, [1009.2450]. 

[20] K. Melnikov and M. Schulze, NLO QCD corrections to top quark pair production in association with 
one hard jet at hadron colliders, Nucl.Phys. B840 (2010) 129-159, [1004 . 3284]. 

[21] G. Corcella et ai, HERWIG 6: An event generator for hadron emission reactions with interfering 
gluons (including supersymmetric processes), JHEP 01 (2001) 010, [hep-ph/0011363]. 

[22] T. Sjostrand, High-energy physics event generation with PYTHIA 5. 7 and JETSET 7.4, 
Comput.Phys.Commun. 82 (1994) 74-90. 

[23] T. Sjostrand, S. Mrenna, and P. Z. Skands, PYTHIA 6.4 Physics and Manual, JHEP 0605 (2006) 026, 

[hep-ph/0603175]. 

[24] G. Bevilacqua, M. Czakon, M. Garzelli, A. van Hameren, Y. Malamos, et ai, NLO QCD calculations 
with HELAC-NLO, Nucl.Phys.Proc.Suppl. 205-206 (2010) 211-217, [1007 .4918]. 

[25] A. Kardos, C. Papadopoulos, and Z. Trocsanyi, Top quark pair production in association with a jet 
with NLO parton showering, Phys. Lett. B211 (201 1) [1101.2672]. 

[26] S. Frixione, Z. Kunszt, and A. Signer, Three-jet cross sections to next-to-leading order, Nucl. Phys. 
B467 (1996) 399^42, [hep-ph/9512328]. 

[27] J. Alwall et ai, MadGraph/MadEvent v4: The New Web Generation, JHEP 09 (2007) 028, 
[0706.2334]. 

[28] K. Hasegawa, S. Moch, and P. Uwer, AutoDipole: Automated generation ofdipole subtraction terms, 
Comput.Phys.Commun. 181 (2010) 1802-1817, [0911 . 4371]. 

[29] T. Binoth et ai, A proposal for a standard interface between Monte Carlo tools and one-loop 
programs, 1001.1307. 

[30] G. van Oldenborgh and J. Vermaseren, New Algorithms for One Loop Integrals, Z.Phys. C46 (1990) 
425^38. 

[31] G. van Oldenborgh, FF: A Package to evaluate one loop Feynman diagrams, Comput.Phys.Commun. 
66(1991) 1-15. 

[32] S. Alioli, P. Nason, C. Oleari, and E. Re, Vector boson plus one jet production in POWHEG, 

1009.5594. 

[33] J. Pumplin et ai. New generation of parton distributions with uncertainties from global QCD 
analysis, JHEP 07 (2002) 012, [hep-ph/0201195]. 

[34] S. D. Ellis and D. E. Soper, Successive combination jet algorithm for hadron collisions, Phys. Rev. 
D48 (1993) 3160-3166, [hep-ph/9305266]. 

[35] S. Catani, S. Dittmaier, M. H. Seymour, and Z. Trocsanyi, The dipole formalism for next-to-leading 
order QCD calculations with massive partons, Nucl. Phys. B627 (2002) 189-265, 

[hep-ph/0201036]. 

[36] J. Alwall et al.,A standard format for Les Houches event files, Comput. Phys. Commun. 176 (2007) 

300-304, [hep-ph/0609017]. 

[37] M. Cacciari, G. P. Salam, and G. Soyez, The anti-k,jet clustering algorithm, JHEP 0804 (2008) 063, 
[0802.1189]. 



-35- 



[38] S. Alekhin, J. Blumlein, S. Klein, and S. Moch, The 3, 4, and 5-flavorNNLO Partonfrom 
Deep-Inelastic-Scattering Data and at Hadron Colliders, Phys.Rev. D81 (2010) 014032, 
[0908.2766]. 

[39] P. Jimenez-Delgado and E. Reya, Dynamical NNLO parton distributions, Phys.Rev. D79 (2009) 
074023, [0810.4274]. 

[40] A. Martin, W. Stirling, R. Thorne, and G. Watt, Parton distributions for the LHC, Eur.Phys.J. C63 

(2009) 189-285, [0901.0002]. 

[41] G. Watt, Parton distribution function dependence of benchmark Standard Model total cross sections 
at the 7 TeV LHC, JHEP m (2011)069, [1106.5788]. 

[42] U. Langenfeld, S. Moch, and R Uwer, Measuring the running top-quark mass, Phys.Rev. D80 (2009) 
054009, [0906.5273]. 

[43] CDF Collaboration, T. Aaltonen et at.. First Measurement of the tt Differential Cross Section 
da/dMff in pp Collisions at^fs= 1.96 TeV, Phys.Rev.Lett. 102 (2009) 222003, [0903 .2850]. 

[44] DO Collaboration, V. Abazov et al. Dependence of the tt production cross section on the transverse 
momentum of the top quark, Phys.Lett. B693 (2010) 515-521, [1001 .1900]. 

[45] DO Collaboration, V. Abazov et al. First measurement of the forward-backward charge asymmetry in 
top quark pair production, Phys.Rev.Lett. 100 (2008) 142002, [0712.0851]. 

[46] CDF Collaboration, T. Aaltonen et al., Forward-Backward Asymmetry in Top Quark Production in pp 
Collisions ats/s= 1.96 TeV, Phys.Rev.Lett. 101 (2008) 202001, [0806.2472]. 

[47] CDF Collaboration, T. Aaltonen et al.. Evidence for a Mass Dependent Forward-Backward 
Asymmetry in Top Quark Pair Production, Phys.Rev D83 (2011) 112003, [1101 . 0034]. 

[48] DO Collaboration, V. M. Abazov et al. , Forward-backward asymmetry in top quark-antiquark 
production, 1107.4995. 

[49] S. Catani and B. Webber, Infrared safe but infinite: Soft gluon divergences inside the physical region, 
JHEP 9710 (1997) 005, [hep-ph/9710333]. 

[50] S. Alioli, K. Hamilton, R Nason, C. Oleari, and E. Re, Jet pair production in POWHEG, JHEP 1104 

(2011)081, [1012.3380]. 

[51] R Z. Skands, Tuning Monte Carlo Generators: The Perugia Tunes, Phys.Rev. D82 (2010) 074018, 
[1005.3457]. 

[52] V. Ahrens, A. Ferroglia, M. Neubert, B. D. Pecjak, and L. L. Yang, The top-pair forward-backward 
asymmetry beyond NLO, 1106.6051. 

[53] J. H. Kuhn and G. Rodrigo, Charge asymmetries of top quarks at hadron colliders revisited, 

1109.6830. 

[54] S. Alioli, J. Fuster, A. Irles, S. Moch, P. Uwer, and M. Vos, A theoretical study of new observables to 
measure the top quark mass at hadron colliders, . 

http : //www . ino . tif r . res . in/MaKaC/ contributionDisplay . py ?contribId=24 Sconf Id=7 9. 

[55] ATLAS Collaboration, Measurement of the charge asymmetry in top quark pair production in pp 
collisions at^fs = 7 TeV using the atlas detector. Tech. Rep. ATLAS-CONF-201 1-106, CERN, 
Geneva, Aug, 201 1. 

[56] CMS Collaboration, Measurement of the charge asymmetry in top quark pair production. Tech. Rep. 
CMS-PAS-TOP-1 1-014, CERN, Geneva, Jul, 201 1. 



-36- 



[57] A. Denner, S. Dittmaier, S. Kallweit, and S. Pozzorini, NLO QCD corrections to WWbb production at 
hadron colliders, Phys.Rev.Lett. 106 (201 1) 052001, [1012 .3975]. 

[58] G. Bevilacqua, M. Czakon, A. van Hameren, C. G. Papadopoulos, and M. Worek, Complete off-shell 
effects in top quark pair hadroproduction with leptonic decay at next-to-leading order, JHEP 02 
(2011)083, [1012.4230]. 

[59] W. Bernreuther, A. Brandenburg, Z. G. Si, and P. Uwer, Top quark spin correlations at hadron 
colliders: Predictions at next-to-leading order qcd, Phys. Rev. Lett. 87 (2001) 242002, 

[hep-ph/0107086]. 

[60] A. Brandenburg, Z. Si, and P. Uwer, QCD corrected spin analyzing power of jets in decays of 
polarized top quarks, Phys.Lett. B539 (2002) 235-241, [hep-ph/0205023]. 

[61] S. Frixione, E. Laenen, P. Motylinski, and B. R. Webber, Angular correlations oflepton pairs from 
vector boson and top quark decays in Monte Carlo simulations, JHEP 04 (2007) 081, 

[hep-ph/0702198]. 

[62] C. OleariandL.Reina, W -\- bb production in POWHEG, JHEP 1108(2011)061, [1105.4488]. 

[63] S. Frixione, E. Laenen, P. Motylinski, B. R. Webber, and C. D. White, Single-top hadroproduction in 

association with a W boson, JHEP 07 (2008) 029, [0805.3067]. 

[64] G. Mahlon and S. J. Parke, Spin Correlation Effects in Top Quark Pair Production at the LHC, 
Phys.Rev. D81 (2010) 074024, [1001 . 3422]. 

[65] P. Uwer, Maximizing the spin correlation of top quark pairs produced at the Large Hadron Collider, 
Phys.Lett. B609 (2005) 271-276, [hep-ph/0412097]. 

[66] DO Collaboration, V. M. Abazov et al. Measurement of spin correlation in tt production using 
dilepton final states, Phys.Lett. B702 (2011) 16-23, [1103 . 1871]. 

[67] K. Melnikov and M. Schulze, Top quark spin correlations at the Tevatron and the LHC, Phys.Lett. 
B700 (2011) 17-20, [1103.2122]. 

[68] DO Collaboration Collaboration, V. M. Abazov et al. Measurement of spin correlation in tt 
production using a matrix element approach, Phys.Rev.Lett. 107 (201 1) 032001, [1104.5194]. 

[69] DO Collaboration Collaboration, V. M. Abazov et al. Evidence for spin correlation in tt production, 
Phys.Rev.Lett. (2011) [1110.4194]. * Temporary entry *. 

[70] M. L. Mangano, M. Moretti, F. Piccinini, and M. Treccani, Matching matrix elements and shower 
evolution for top-quark production in hadronic collisions, JHEP 01 (2007) O13,[hep-ph/0611129]. 

[71] ATLAS Collaboration, Measurement of spin correlation in tt production from pp collisions at 
^/s = l TeV using the atlas detector. Tech. Rep. ATLAS-CONF-2011-117, CERN, Geneva, Aug, 
2011. 

[72] W. Bernreuther and Z.-G. Si, Distributions and correlations for top quark pair production and decay 
at the Tevatron and LHC, Nucl.Phys. B837 (2010) 90-121, [1003 .3926]. 

[73] S. Ahoh, K. Hamilton, and E. Re, Practical improvements and merging of POWHEG simulations for 
vector boson production, 1108.0909. 



-37- 



